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ABSTRACT 

We review the effective Lagrangian of the Higgs penguin in the Standard Model and its minimal 
supersymmetric extension (MSSM). As a master application of the Higgs penguin, we discuss in 
some detail the B-meson decays into a lepton-antilepton pair. Furthermore, we explain how this 
can probe the Higgs sector of the MSSM provided that some of these decays are seen at Tevatron 
Run H and B-factories. Finally, we present a complete list of observables where the Higgs penguin 
could be strongly involved. 
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1 Prologue 



Glashow and Weinberg in their seminal paper P pointed out that Flavour Changing Neu- 
tral Currents (FCNC) are suppressed "naturally" if all the down-type quarks acquire their 
masses through their coupling to the same Higgs boson doublet, say Hd, and all the up- 
type quarks through their coupling to a second Higgs boson doublet, say if^. These two 
doublets may [Standard Model (SM) case] or may not [Two Higgs doublet Model type-II 
(2HDM) case] be charge conjugates of each other. In general, if no symmetry considera- 
tions are assumed at tree level, a departure from the above rule leads to severe enhancement 
of K — K and/or B — B mixing not seen at the experimental data |2j. In the Minimal 
Supersymmetric Standard Model (MSSM) the Glashow- Weinberg scheme is naturally re- 
alized due to the holomorphicity of the superpotential; Higgs FCNC processes appear only 
at loop level when the holomorphicity is violated by finite radiative threshold corrections 
due to the soft SUSY breaking interactions In this brief review, we shall present the 
effective Lagrangian of the Higgs boson FCNC's in the SM and MSSM, which we term 
"Higgs penguins", and discuss their significance in B-, K- meson and r-lepton physics. 



2 The Higgs penguin 

The term "Higgs penguin" is used here in analogy to the well known Z-penguin, to denote 
H — f — f loop induced flavour transitions, where /, /' are either quarks or leptons. For B- 
or K-physics experiments, the relevant Higgs penguin is the one with down quarks (/ = d) 
in the external legs, drawn schematically in FigCl 



/ /' 




Figure 1: The Higgs penguin. The blob consists of loop corrections due to Standard or 
Beyond the Standard Model particles. 
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2.1 The Standard Model Higgs penguin 



The flavour changing Higgs vertex H — d — d! m the SM was first calculated in . It 
came rather as a surprise to observe that the coupling of the resulting Higgs H — d — d' 
penguin vertex is proportional to md(d')/Mw instead of being suppressed by {md{d') / Mw)^ ■ 
Denoting the up and down 3x3 diagonal quark mass matrices with Mu and Md? and with 
V the Cabbibo-Kobayashi-Maskawa (CKM) matrix, the SM one-loop effective Lagrangian 
for quark flavour changing interactions with the Higgs boson reads 

= " 2^ ^ ^ ( + ^'^'^^ ) ^' ' ^^-^^ 

where Pl{r) = [1 - (+)75]/2, and 

L_ _ 9w yt ^^^u y R_ _ / L _ Nf 2) 

^Hdd'{d^d')- ^ (I67r2) ' - \^H,dd' ) ■ \^-^> 

The tree level flavour diagonal quark Higgs couplings are g^j^ = {^^.gd)^ ~ ^- -^^^ 
this calculation the masses and the momenta of the external particles have been set to 
zero. This simple result arises from the sum of 10 diagrams in gauge (6 Higgs pen- 
guins and 4 self energy diagrams) shown in Fig.Q. In a general i?^ gauge, the dia- 
grams (c),(e),(f),(g),(h),(j),(k) are divergent. The Glashow-Iliopoulos-Maiani (GIM) mech- 
anism IHj removes divergences from diagrams (e) and (f) aswell as from (g) and (j). The 
remaining divergences in diagrams (c), (h) and (k) cancel each other and the finite sum, 
in the case of a constant background Higgs field^, is given in Eqs. ()2.1l2.2p . Although 
it seems straightforward to calculate the 10 diagrams in Fig.®, it took many years for 
gauge dependencies and renormalization scheme ambiguities to be clarified in the liter- 
ature. The result of Refs. [HE] given in Eqs. ()2.1l2.2|) was finally confirmed by several 
groups [HIHlinilinilllllllllllll and different methods. 

^For non zero Higgs boson mass, the coupling g^V t,,, . , is in general written as: 



''Hdd'{d^d' 

„ , . _ 3w 

^Hdd'(d^d') ~~16^ 



L _ 5„vt/(x,y)V, (2.3) 



with X = -fjT- and y — -j-^-, and IS] 



(~4(T^ + 2(1^ - 8(1^ + 8(1^) ■ ^'-^^ 

Numerically /(a; 0, j/) = 0, /(4,0) = 3, /(4,2) = 3.7, /(4,4) = 4.4 and /(4, 100) = 38.5. Thus, 
for a Higgs mass around the electroweak scale, the approximation made in deriving Ea. (|2.2|l is good. 
However, as the Higgs mass goes over the TeV scale, the ^^g^n^^^^,-^ enters into the non-perturbative 
regime. Furthermore, for applications it is the ratio J{x,y)/y which appears in the physical amplitudes. For 
example, for a Higgs mass Mh = 2.6 GeV ^y = 0.001 ^ /(4, y)ly = 3000, and limy^oo /(4, y)ly = 0.36. 
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Figure 2: SM contributions to the Higgs penguin. The label u denotes all the up-type quarks 
{u,c,t), and the charged Goldstone modes. 



For B^- or i?s-nieson initial state [d = b and d' = d or s), the dominant contribution 
in Eq. ()2.2|) comes from the top quark in the loop while for Kaon initial state {d = s and 
d' = d) it comes from the charm quark in the loop. Comparing the Higgs penguin prefactor 
of Eq. ()2.H) and Eq. ()2.2|l to the Z-boson penguin [d Zy, d') one in Ref. JH], we obtain that 
the former is about rrid/Myy times smaller than the latter, with being the mass of the 
heaviest external quark. 
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2.2 The Supersymmetric Higgs penguin 

For the construction of the supersymmetric version of the SM Higgs penguin two approaches 
have been devised: the Feynman diagrammatic approach and the Effective Lagrangian 
approach. Both have advantages and disadvantages and it is worth describing briefly both 
approaches here. 

2.2.1 Feynman Diagrammatic Approach 

This approach is a straightforward (but somehow tedious) calculation of the diagrams shown 

in Fig.Q. It needs to be supplemented with the so called tan /3-resummation procedure. 

The main advantage is that within this approach effects beyond the SU{2) x f/(l) limit 

are automatically included. This approach has been followed by various groups (usually 

through a specific application) [TTllTRllTTIll^lTni^ l^ . The 17 diagrams in Fig.® are just 

supersymmetrizations of the SM diagrams in Fig.©. To our knowledge they have not yet 

been calculated in full detail ^. The latest calculation ^H] is restricted only to results that 

have been derived in the large tan/3 approximation with the low tan/3 regime remaining 

unexplored although one does not expect dramatic effects for heavy (TeV) supersymmetric 

particles^. Some remarks on the diagrams in Fig.Q are in order here: i) the bulk of the 

corrections to the Higgs penguin arise just from the finite part of the self energies (g,j) 

with charginos and up-squarks in the loop; they scale like tan^ (3. They do not vanish 

even (and especially) in the case of degenerate squark masses, and the MSSM corrections 

to the Higgs penguin do not decouple even if the squark masses are at multi-TeV scale. 

Diagrams (d),(e),(f),(i),(g) and (j) with the charged Higgs boson (two Higgs doublet model 

contributions [211123^1 ) are also non zero but are sub dominant since they scale like 

tan/3. Ultraviolet Infinities of the diagrams (g,j) cancel against the vertex diagram (b). ii) 

Diagrams (a),(b),(c),(h) and (k) with neutralinos or gluinos and squarks in the loop arise 

only for non degenerate squark spectrum due to renormalization group running effects or 

Planck scale squark flavour changing insertions. In this case, the gluino contributions can 

compete with the chargino ones and even cancel each other (this will be clear in a minute). 

The neutralino diagrams in (a,b,h,k) play rather a subdominant role except for the case 

of the exact cancellation of the gluino and chargino contributions [201 • Within the 

After posting the paper, the author was informed that in the Umit of vanishing external momenta of 
quarks and the Higgs particle, a complete calculation i.e., no expansion in tan/3, exists in a numerical code 
in Ref. jTH]. 

"^Non-tan /3 enhanced electroweak corrections tend to decouple rapidly as the supersymmetric particles 
approach the TeV scale. 
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Figure 3: General and complete MSSM contributions to the Higgs penguin. The labels u, d 
denote all flavours for up and down type squarks with a=l,2 being their mass eigenstates. 
Charginos are labeled with x~, « = 1, 2, neutralinos with i = 1...4 and the gluino with g. 
is the charged Higgs boson. 

Feynman diagrammatic approach, tan /3-resummation effects have been added in Ref. 
for all diagrams in Fig.Q apart from the gluino and neutralino ones. In addition, SUSY 
CP- violating effects were taken into account in j22j . 
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2.2.2 Effective Lagrangian Approach 



This approach allows for a simpler derivation of the dominant corrections to the Higgs 
penguin (including resummation of tan (3) in the large tan /? ^ 40 regime but can not easily 
implement SU{2) x U{1) electroweak symmetry breaking effects which are relevant for 
small values of tan/3. This approach has been followed in Refs. [^17711^1^1^011^ 1^. 
The general resummed Higgs penguin, which includes CP-violating effects in the MSSM, 
is given by PU] : 

3 

^H.M' = -^J2H^d{M^ z'urM Pl + zIm' Md'Pii ) rf' , (2.5) 
in the notation of Eq. fj2.1jl . with 



On 02i ^ iOsi 



g2,j.'=( )^ (2.6) 



cos P sin P sin P cos /? 

where R is the 3x3 dimensional mass matrix which resums all the tan /3 enhanced finite 
threshold effects. It is given by : 

R = 1 + Egtan/5 + Eutan/5|hu|^ + ... , (2.7) 



where hu is the diagonal up- Yukawa couplings, and up to small corrections is : hu = ■ 
Furthermore, Oij is a 3 x 3 orthogonal matrix, which transforms the Higgs boson fields 
from their weak to their mass eigenstates and accounts for the CP-violating Higgs mixing 
effects jnS]- In addition. Eg and Eu are finite threshold effects induced by the diagrams 
in Fig.(jD). The ellipses in ()2.7p denote additional (generically sub-dominant) threshold 
effects like, for example, the bino-wino contribution which corresponds to the neutralino 
contribution diagram in Fig.®. Eg and Eu are in general complex 3x3 matrices but, for 
simplicity, here we will assume that they are diagonal (the reader may consult Ref. ^2Qj for 
more details at this point) : 

Eg = l|im;;.*/(m|^,m|^,Kr)^-^^-l^e-(^''-^^«), (2.8) 



The limit in Eqs. ()2.8l2.9|) is at degenerate supersymmetric masses, Msusy = = |/^| = 
mj^ = mj^ = \Au\ = mnj^ = m^^ and (j)^, (pg, (pAu the phases of the higgsino mixing 
parameter /i, the gluino mass rrig, and the trilinear supersymmetry breaking coupling Au, 
respectively. The loop integral I{x, y, z) is given in [nHj. Even in this limit, Msusy ^ My/ 
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Figure 4: Non-holomorphic radiative effects on the down-quark Yukawa couplings of the 
form [hd (ii$2*'^fi] induced by (a) gluinos and (b) higgsinos. The Higgs field $2 ^■^ the one 
which couples to the up-type quarks in the superpotential, i.e., [huW^^^^^ij]. From Ref. 



a nightmare scenario for the LHC, the Higgs penguin amphtude remains significant and 
does not decouple; it depends on the various CP-violating phases and details of the Higgs 
sector [mixing matrix {Oij)] and the experimentally nearly-known quark sector. 

Equations ()2.5I2.6I2.7|) represent the Higgs penguin amphtude in the MSSM. In order 
to qualitatively understand the behaviour of the MSSM Higgs penguin, it is instructive to 
take the limit of large tan/3, i.e., 1/ cos/3 ~ tan/3, sin/3 ~ 1 and apply the limits ()2.8p2.9|l 
in Eqs. fl2.5f2.fif2.7|l . Then, we obtain: 



On + i Osi] e-'^^^+'t'^u) tan2 /? 



^+<^9) tan/3 



Q-i{<i>^+(t>g) _)_ 



9I M2 
647r2 ^ 



E/)) tan/3 



^Hidd' )^ 



V 



(2.10) 



By looking at Eq. ()2.1()j) . we conclude the following : 



The MSSM Higgs penguin is enhanced by two powers of tan/3 with respect to the 
SM one. Compare Eqs ^(TY^ and 

The MSSM Higgs penguin is in general a complex number due to the additional 
super symmetric CP- violating phases, ipAu-, 4>ti- In the CP-invariant limit the coupling 
dd' {d^d') either pure real or pure imaginary. 

The MSSM Higgs penguin's largest correction, Eq. ()2.9|) . depends on the soft super- 
symmetry breaking trilinear parameter, Au- Therefore, it is enhanced only in the 
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case of the maximal stop- mixing case, At ^ Mw, which is favoured also by MSSM 
Higgs searches. 

• The limit [1 + Eg tan j3] in the master formula Eq. ()2.6|) . or equivalent the approx- 
imation [1 + ^ Q-A't>i^+<f>g) tan/5] — i> in Eq. ()2.10|) . is attainable and not singular. In 
this limit the Higgs penguin diagram goes to a constant value (and not to infinity) due 
to the GIM mechanism. This is a result of the general resummed effective Lagrangian, 
equipped in jSU], that integrates and improves earlier constructions [2111211123121] • 

• The MSSM Higgs penguin will dominate over all the box- or penguin-type supersym- 
metric contributions to the physical observables since the latter vanish due to the 
super GIM mechanism in the limit of the degenerate squark masses pHl35j . 

• The resummation matrix R controls the strength of the Higgs-mediated FCNC ef- 
fects. For instance, if R is proportional to unity, then a kind of a GIM-cancellation 
mechanism becomes operative and the Higgs-boson contributions to all FCNC ob- 
servables vanish identically. Furthermore, not only the top quark, but also the other 
up-type quarks can give significant contributions to FCNC transition amplitudes, 
which are naturally included in ()2.5|) through the resummation matrix R. All these 
effects are computed explicitly within physical observables in Ref. |3D|. 

• The Lagrangian in Eq. ()2.5|) captures the bulk of the full corrections of Fig. Q at large 
tan /5 ^ 40 and is limited for supersymmetric breaking masses Msusy much heavier 
than the electroweak scale. Additional electroweak corrections (after the SU{2) x U (1) 
gauge symmetry breaking) have been calculated in [HI]. 

We will now turn to the physical applications of the Higgs penguin. 

3 Applications 

There is a vast amount of applications one can think of. Basically, the Higgs penguin 
participates in those physical observables where the Z-penguin does so. In this section we 
shall discuss, in a rather qualitative way, the effect of the Higgs penguin [H — b — d{s)] 
on B-, K-meson and r-lepton physics observables. The models under consideration will 
be the Standard model and its minimal supersymmetric extension, MSSM, with R-parity 
symmetry. 
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i3(Channel) 


Expt. 


Bound (90% CL) 


SM prediction 


Bs e+e" 


L3 ng 


< 5.4 X 10"^ 


(8.9 ±2.3) X 10-1^ 


Bs fi^fj," 


CDF gni 


< 9.5 X 10"^ 


(3.8 ±1.0) X 10-9 


B, ^ r+r- 


LEP isni 


< 0.05 


(8.2 ±2.1) X 10"^ 


Bd ^ e+e- 


Belle ISH 


< 1.9 X lO"'^ 


(2.4 ± 1.4) X 10"^^ 


Bd -> 


Belle HH 


< 1.6 X 10-^ 


(1.0 ±0.6) X 10-1° 


Bd ^ r+r- 


LEP ISHj 


< 0.015 


(2.1 ± 1.2) X 10-^ 



Table 1: The experimental status and the SM predictions for the branching ratios B(Bq 

). The error in the Bg branching ratios mainly originates from the uncertainty in 
fBs = 230 ± 30 MeV and the error in the B^ branching ratios corresponds to the 

uncertainties in fs^ = 200 ± 30 MeV f^J and \Vtd\ = 0.040 ± 0.002 linearly added. This is 
the updated version of the Table 1 presented in Ref. 

3.1 The Master Application : Bq l^l^ and implications 

It is rather compulsory to discuss first the i?q-meson [q = d, s) decay to a lepton (/ = e, fi, r) 
pair. The decay Bg is mediated by the Higgs penguin, the Z-penguin and box 

diagrams [12] . It has been extensively discussed in the literature, both in the SM and in 
the MSSM jnil case (for other models see Ref. jUj). None of the processes Bq — > have 
been seen so far, and the experimental bounds on these "very rare" B-decays are listed in 
Table 1. The decays Bq —>■ are interesting probes of physics beyond the SM since the 
best bound we have so far, that on B{Bs fi'^fi^), is three orders of magnitude away from 
the SM expectation 40]. In addition, the decay Bs^d —* fi'^ ^~ is one of the experimentally 
favoured and has a bright future : Tevatron Run II with luminosity of 2 fb-^ has a single 
event sensitivity of B{Bs — > = 1.0 x 10^^ (the background from Run I is roughly 

one event): 10 events at CDF with C = 2 fb^^ means a B{Bs — > yU^/i") — 10^''. A recent 
analysis |HHHT] showed that CDF can discover Bg fJ'^f^~ in Run lib with an integrated 
luminosity of 15 fb"^ if B{Bs /U+/i-) > 10"^ If B{Bs fi+f^-) turns out to be SM-like 
(see Table 1), then only LHC will be able to measure it |55|. 

Motivated by the above discussion, we are now going to see why and how one can 
achieve branching ratios for the decays Bg —>■ l^l^ that could appear soon at the Tevatron 
and B-factories. To this end, we shall follow the effective Lagrangian technique employed 
in Ref. [SO]- Neglecting contributions proportional to the lighter quark masses nid^s, the 
relevant effective Hamiltonian for AB = 1 FCNC transitions, such as 6 ^ qi^i^ with 
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q = d, s, is given by 



n. 



AB=1 

eS 



where 



AGf 



CsOs + CpOp + CoO 
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(3.11) 



Of 



o 
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167r2 
167r2 
167r2 



(3.12) 



By making use of our master resummed Higgs penguin effective Lagrangian Eq. p.5|2.6|2.7j) . 
the Wilson coefficients Cs and Cp (in the region of large values of tan/3) read: 



Cs 



E 



Ml 



2'irmi 1 



a, 



em ^ife^4 i=l 



9H,qb9H,l 

Ml 



(3.13) 



with Cio = ~ sin^6' ( i67 GeW ~ '^''^ E^l ^ud mt{mt) the running MS* top quark mass. 

denotes the neutral Higgs boson masses, Mh^ < ^ Mh-^. Cio is actually the 
dominant contribution in the SM case and originates from the Z-boson penguin and the 
box diagrams [T3]. In addition, the reduced scalar and pseudoscalar Higgs couplings to 
charged leptons (7^^^ in Eqs. fj3.13|) are given by 



cos j3 



p 

9Hd 



tan P O. 



3i 



(3.14) 



where loop vertex effects on the leptonic sector have been omitted as being negligibly small. 
The reader must have already noticed that, at the large tan/3 regime, we have g'l^^ — 
On X tan/3 and thus from Eqs. ()2.1()|3.13p3.I^ we obtain that the scalar and pseudoscalar 
Wilson coefficients Cs,p grow like tan'^ (3. The branching ratio for the meson decay to 
acquires the simple form [T^ l^ 



i3(5j 



£+r) = 

^^"^ Mb tb WftVr^ 



'tq\ 



'1- 



4m„ 



M 



Ba 



1 - 



Ami 



I P9|2 



+ 



(3.15) 



+ 



2m^Fj|2 
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Figure 5: The B{Bs — >■ ^'^fi~) vs. the Heaviest Higgs boson mass in the MSSM for 
tan P = 50 and rrit = 175 GeV. The shaded area shows the SM prediction for B{Bs /^^/U^) 
(see Table 1). The excluded area from Tevatron CDF \49\l is also displayed. 



where r^,, Mb^ is the total hfetime and the mass of the Bq meson, respectively, and 

- - I "y^- ^ • C.. . (3.16) 

Fs^p is roughly (m^/M^) tan'^ /5 times bigger than F^ and thus dominate for large values 
of tan /3. Other electroweak SUSY box and Z-penguin contributions grow with at most 
two powers of tan j3 and thus are subdominant in this region. Numerical values for the 
parameters entering Eq. ()3.15fH.l(ij) can be found in PDG [21 and in the caption of Table (1). 
The Higgs matrix elements [Oij] in the MSSM can be obtained from the numerical code 
CPsuperH [52] • In general, we have OnjO^i ^ 1, and thus from Eqs. ()3.13|) we conclude 
that the ratio B{Bq depends only on the heaviest Higgs boson mass, Mh^- In the 

simple case of equal squark masses, and by employing Eqs . ()2 . 1 ( ) p3 . 1 3p3 . 1 4|3 . 1 5|3 . . the 
reader can convince himself that the experimental bounds in Table 1 are attainable in the 
MSSM. In general cases of squark non-degeneracy one should use the master formula given 
in Eqs.((23ESD. 
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3.1.1 Implications: probing the Higgs sector with Bg^d l^'^l^ at Tevatron 
and B-factories 

As we show, in the MSSM B{Bs /W^At ) is enhanced by six powers of tan (3 and suppressed 
by four powers of the heaviest Higgs boson mass, Mn-i- On the other hand, tan/5 has 
a theoretical upper bound (around 50) coming from the perturbativity of the Yukawa 
couplings up to the Grand Unification (GUT) scale. Therefore, we can use a future possible 
evidence for Bg — > fi^fi^ to set an upper bound on M^^ ■ Having all the machinery at hand 
from the previous section, we do this exercise here. In order to illustrate our argument, let 
us envision the following situation: 

• Tevatron Run 11 finds Bg fJ'^fJ'~ 

• BaBar and Belle (or a high luminosity upgraded B-factory) confirm the equation 
B{Bd — \^\'^B{Bs — > /i"'"/i~), and no indication is found for non- minimal 
flavour structure 

What do these hypothetical facts imply? We can use a high statistics scan over the MSSM 
parameter space in the region where the soft breaking masses are below 2.5 TeV and 
the trilinear couplings below ^ 5 TeV. We assume that the Lightest Supersymmetric 
particle (LSP) is stable and require it to be neutral. No new flavour structure other than 
the CKM is assumed. Under these assumptions, we plot in Fig.(jHI) B{Bs f^^f^") vs. 
Mhs- The envelope contour is well approximated by 

, , 7 /tan/3\6 /550 GeV\4 „ 

B(B. ^ = 5 X 10- (^) (^^) + 8 X 10-» . (3.17) 

If CDF Run Ha sees 50 events for the Bg ^~ [that means B{Bg — > ~ 5 x 10"''] 

then Mh^ will be less than 550 GeV for all tan/5 values less than 50. The fit-equation 
(K117|l therefore sets an upper bound on Mh-j,, a very useful result indeed. We should note 
here that various phases, or non-minimal flavour structure in the squark sector, can alter 
the above result. A more detailed analysis is beyond the scope of this review. 

3.2 Other Applications of the Higgs Penguin 

The Higgs penguin applications can be classified into the six categories depicted in Fig.®: 

• Quark single Higgs penguin decaying into a pair of leptons (1) or quarks (q) of the 
same family [FiglHfA)]. This category includes the master application B —>■ 1 1 [HZ! 
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Figure 6: Classification of all the applications of the Higgs penguins. The vertex-dot indi- 
cates that the coupling is proportional to tan (3 for the weak isospin (—1/2) quarks q = d,s,b 
and leptons I = e, fi, r. It is suppressed by cot (3 otherwise. The down quark Higgs penguin 
(blob with d,d' quarks) grows with two powers o/tan/5 in the MSSM. 



that we have already discussed in the previous section. It leads also to contributions 
to the semileptonic B-decays B X II jl9 |l!MIF7j and K II. Forward-backward 
asymmetries m. B ^ Xrr, B ^ nil and K it 1 1 influenced by the Higgs penguin 
have been studied in |2H], [SHI and [Hni, respectively, and the lepton polarization 
asymmetry m. B ^ II oi K ^ II decays has been studied in [;61j and jH^! respectively. 
The decay into a quark pair d' — ^ dqq leads to a number of processes like K — > tttt jHOj . 
B ^ nn,B ^ (pKs, B J/^ K,B ^ Ktt etc. In the MSSM, the amplitude grows 
like tan^ /3 if g = ci, s (or / = e, /x, r), or with tan/5 ii q = u, c. We have seen already 
that the Higgs penguin is in general complex [see for example, Eqs. ()2.1|2.5p ]. As 
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a consequence, one can construct important CP-observables like, for example, the 
leptonic CP-asymmetry [HSll^ lHUj. 

Lepton Flavour Violating (LFV) single Higgs penguin decaying into pair of leptons 
(1) or quarks (q) [FigE^B)]. The LFV Higgs penguin has been calculated in the SM 
in Ref. [Hi]. In the MSSM, the method for calculating the Higgs penguin, presented 
in Ref. jHOI and summarized here, may be easily extended to consistently account 
for charged-lepton flavour violating effective Lagrangian. The r-decays into three 
leptons r ^ (e, /i) H [HSj, or r-decays into a lepton and a pair of quarks r (e, ij)qq, 
like r ^r], belong to this category. The amplitude of these decays are neutrino 
mass generation mechanism dependent. For example, in the see-saw MSSM, t ^ 31 
decays [OniinZllSSl are enhanced by three powers of tan/5 with the same happening 
in the r — >• Iqq decays or ftq eq amplitude (/i — e conversion in the nuclei) 
recently studied in [70]. On the contrary, decays with r — > luu grow with only one 
power of tan /5. 

Quark-Lepton Flavour violating double Higgs penguin [FigEfC)]. This is nothing else 
than the combination of the cases (A) and (B). The processes B —>■ W jHIj, B XW 
and T (e, ^) K, etc belong to this category. They start formally at the two- 
loop level and they grow with four powers of tan (3. It is very interesting to notice 
that the double penguin (C) leads to invisible (at the detector) B-meson decays if 
the final particle state is a pair of neutrinos (u) of different flavour. This would be a 
unique signature at B-factories; we believe that it should be further studied from the 
theoretical point of view and searched for experimentally. 

Quark double Higgs penguin [FigE^D)]. Such processes lead to contributions to 
B — B OT K — K mixing, and in the MSSM dominate over all other contribu- 
tions for large tan/3. Although there is no similar analysis in the SM, following 
our discussion above, one expects that the diagram (D) is suppressed by a factor 
{rrib/Mw)'^- In the MSSM, the dominant contribution to the Bs^d — Bg^d mixing is 
proportional to mbfns^ta.n'^ P / M]^^ [7^l3Up31j . and to the K — K mixing proportional 
to mdnis^d^ai]^ [3 / M'jj^ jSU]. In addition, a complex double penguin leads to correc- 
tions on the CKM angle-/? originated from soft supersymmetry breaking complex 
parameters [HH] . 

Quark-Gluon double Higgs penguin [FiglHl^E)]. This is another class which is a com- 
bination of two penguins: the triangle Higgs-Gluon one (effective Lagrangian in the 
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SM and in the MSSM can be found in [7S]) and the Quark Higgs penguin reviewed 
here. They contribute to decays B {X)gg and they scale at most hke three powers 
of tan/3. No study of this Higgs mediated process exists so far in the hterature. 

• Quark-Photon double Higgs penguin [FiglHl^F)]. A combination with the photon tri- 
angle Higgs vertex - the main Higgs decay mode search at the LHC jZH]- It leads to 
processes like B —>■ (X)77 and scales at most like tan^ j3 when bottom quarks run in 
the photon triangle. No study for this amplitude exists in the literature. 

All the above applications should be correlated with each other if they dominate in the 
corresponding processes. One example is the correlation between B{Bs /U^/W ) and 
Bg — Bg mass difference, AMg, found in [72]. This correlation is valid in a specific scenario 
of squark mixing and is lost in a general case pm i31j. Correlations among the full set of 
observables in FiglHl would be an interesting probe of the MSSM squark flavour structure. 
As we can see, there is still a lot of work to be done along that direction. 

4 Epilogue 

In this brief review, we discussed the effective Lagrangian of the Higgs penguin both in 
the SM and in the MSSM, which are presented in Eqs. ()2.H) and ()2.5p . respectively. The 
SM Higgs penguin exhibits a non-decoupling behaviour proportional to V"'' (M^/M^) V, 
but it is suppressed relative to the Z-penguin by a factor of ~ rrid/Mw- The MSSM Higgs 
penguin exhibits also a non-decoupling behavior but is further enhanced by two powers of 
tan /3, and dominates over the Z-penguin if tan (3 is large. We have in addition presented 
the applications of the Higgs penguin in B-meson, K-meson and r-lepton physics. The 
master application is the decay Bg —>■ . In the SM, the Higgs penguin contributes 
negligibly to B{Bq if the LEP experimental bound for the Higgs boson mass is 

taken into account. On the other hand in the MSSM, it is enhanced by three orders of 
magnitude mainly due to the tan/3 enhancement shown in Eq. ()2.10j) . The B-decay mode 
Bs4 — > fi~^fi~ is currently under high search priority at Tevatron and B-factories. As we 
explained in section 3.1.1 one of the most interesting aspects is the fact that a possible 
evidence for Bg fi^fi^ can probe the Higgs sector of the MSSM. Finally, we unfolded an 
exhaustive list of applications in Section 3.2, some of which are novel and not thoroughly 
studied even in the SM. We strongly believe that the Higgs penguin amplitudes deserve 
further theoretical and experimental investigation as their processes may soon leave their 
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footprints at Hadron colliders and B-factories. 
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